Generalizing the generalized Chaplygin gas 
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The generalized Chaplygin gas is characterized by the equation of state p = —A/p a , with a > — 1 
and w > — 1. We generalize this model to allow for the cases where a < —1 or w < — 1. This 
generalization leads to three new versions of the generalized Chaplygin gas: an early phantom 
model in which w <C — 1 at early times and asymptotically approaches w = — 1 at late times, a late 
phantom model with w ~ — 1 at early times and w — > — oo at late times, and a transient model with 
w ~ —1 at early times and w — > at late times. We consider these three cases as models for dark 
energy alone and examine constraints from type la supernovae and from the subhorizon growth of 
density perturbations. The transient Chaplygin gas model provides a possible mechanism to allow 
for a currently accelerating universe without a future horizon, while some of the early phantom 
models produce w < — 1 without either past or future singularities. 



I. INTRODUCTION 

The universe appears to consist of approximately 30% 
dark matter, which clusters and drives the formation of 
large-scale structure in the universe, and 70% dark en- 
ergy, which drives the late-time acceleration of the uni- 
verse (see Ref. for a recent review, and references 
therein) . Numerous models for the dark energy have been 
proposed; in one class of models, the dark energy is sim- 
ply taken to be barotropic fluid, in which the pressure p 
and energy density p are related by 



P = f(p)- 



(1) 



One of the first cases to be investigated in detail was the 
equation of state 



P = wp, 



(2) 



where w is a constant 0, Q . For an equation of state of 
this form, the density of the dark energy depends on the 
scale factor, a, as 



p oc a 



-3(^+1) 



(3) 



A more complex equation of state arises in the Chap- 
lygin gas model 0,[H, for which the equation of state has 
the form 



A 

1 

P 



(4) 



where A is a constant, leading to a density evolution of 
the form 



(5) 



In this model, the density interpolates between dark- 
mattcr-likc evolution at early times and dark-energy-like 
(i.e., constant density) evolution at late times. Thus, the 
Chaplygin gas model can serve as a unified model of dark 
matter and dark energy. (The dark matter sector of this 
model may have problems with structure formation 



although see the discussion in Refs. [1,13 for an alternate 
viewpoint). 

The Chaplygin gas model was generalized by Bento, et 
al., E3 to an equation of state of the form 



Pgcg 



A 



gcg 



(6) 



where A and a are constants. This equation of state 
leads to the density evolution 



Pgcg 



A 



B 



,3(1+0) 



l/(l+a) 



(7) 



Again, the density evolution in the generalized Chaply- 
gin gas models changes from p oc a~ 3 at early times to 
p = constant at late times. Note that such an equation 
of state can also be modelled as a dissipative matter fluid 
where the dissipative pressure is proportional to the en- 
ergy density; a number of exact solutions for this model 
have been discussed in Refs. [TlLIT^|. 

In equation one normally takes a > — 1, while A 
is chosen to be sufficiently small that w > — 1 for the 
generalized Chaplygin gas. In this paper, we relax these 
constraints and examine Chaplygin gas models described 
by equation © for which a and w can lie outside this 
range. In this case, the Chaplygin gas no longer serves 
as a unified model for dark matter and dark energy, but 
it can be taken to be a model for dark energy alone. 
This sort of generalization for the special case of a = 
1 was previously undertaken by Khalatnikov 13]. The 
models presented here can be derived as special cases 
of the more generic models examined in Refs. 

El El, 

and can also arise in the context of fc-essence models 03 . 
These possible generalizations are mentioned in Ref. [TtJ , 
although only the case a > — 1, w < — lis explored in 
detail. In addition to noting some interesting features of 
these models, our main new result is a set of observational 
constraints on these models. 
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FIG. 1: The energy density p of the generalized Chaplygin 
gas (normalized to the present density po) as a function of 
the scale factor a (taken to be 1 at present). Solid curve is 
Case 2(a) (early phantom with a > 0: A s — 1.1, a = 0.3), 
dashed curve is Case 2(b) (early phantom with — 1 < a < 0: 
A s — 1.2, a — —0.95), dashed-dot curve is Case 3 (transient 
GCG, with As — 0.9, a = —1.5) and dotted curve is Case 4 
(late phantom, with A s — 1.1, a = —1.1). 



II. GENERALIZED CHAPLYGIN GAS AS A 
DARK ENERGY COMPONENT 

The equation of state for the generalized Chaplygin 
gas (GCG) is given by equation In order to produce 
a negative pressure and give the currently observed ac- 
celeration of the universe, equation © must have A > 0, 
so we will confine our attention to this case. We assume 
further that a ^ —1, since for a = —1, the GCG is 
equivalent to the dark energy fluid described by equation 

©■ 

Note that a particular choice of A and a does not 
uniquely determine the GCG model; one also needs to 
specify, for example, the value of p at a particular red- 
shift. Integration of the energy conservation equation 



T£ = for a ^ — 1 results in 



Pgcg = Pg Cg0 [A s + (1 - A,)(l + Z)^ i+a >] 



3(l+a)il/(l+a) 



where 



A s = A/p 1+a 



gcgO' 



(8) 



(9) 



with Pgcgo being the present value of p gcg - The choice of 
A s and a does uniquely specify the GCG model. 

Since we are considering the GCG as a dark energy 
candidate, the expression for the Hubble parameter near 
the present becomes 



FIG. 2: As Fig. 1, for the generalized Chaplygin gas equation 
of state parameter w, as a function of the scale factor a. 



H- 



^mo(l + z ) 3 + Qgcgo[A s 



+ (1- A s )(l + ^) 3(1+Q) ] 1/(1+Q) 



(10) 



where we have assumed a normal dust-like dark mat- 
ter component with present density parameter f2 m o, and 
QgcgO is the present value of the density parameter for 
the GCG component. We have neglected the radiation, 
which makes a negligible contribution to the total den- 
sity around the present epoch. We assume a flat universe: 
Qgcgo = 1 — fi m 0' The equation of state for the GCG fluid 
is 



A s + (l-A s )(l + z)^+«)' 



Taking z — in this equation, it is clear that 



wo 



-A* 



(11) 



(12) 



where wq is the present-day value of w for the Chaplygin 
gas. Equation (|12f) gives the physical significance of A s . 

The behavior of the GCG can vary significantly, de- 
pending on the values of A s and a. Our aim in this 
paper is to explore all of these possible behaviours. 

We consider first two trivial cases. For A s = 1, the 
GCG behaves exactly as a cosmological constant at all 
times for all values of a. For A s = 0, the GCG behaves 
as a standard pressureless (w = 0) dust component for 
at all times for all values of a. 

Now consider the four non-trivial cases corresponding 
to the pair of choices A s < 1, A s > 1, and a < — 1, 
a > -1. 



Case 1) < A s < l,a > -1 {Standard GCG). At 
previously discussed, in this region the GCG behaves as 
a pressureless dust component at early times, evolving 
asymptotically to a de-Sitter regime (w — » —1) at latf 
times. Hence the GCG in this parameter region car 
act as a unified model for dark energy and dark mat- 
ter (UDM). This is the standard model of the general- 
ized Chapl ygin gas, p reviously explored in great detai 

[TSIIllillailia. 

The remaining three cases correspond to "new" ver- 
sions of the generalized Chaplygin gas. 

Case 2) A s > 1, a > -1 (Early Phantom GCG). Ir 

this case, the GCG acts as a phantom component wit! 

w < — 1 at all times, but w asymptotically approaches — ] 

at late times. Hence, one has an early phantom behavioui 

in this region of parameter space. In this case, p gcg — C 

i 

at 1 + zi, = [A S /[A S — l)] 3l - 1 + a ) and then p gcg grows tc 
become a constant at late times. 

The behavior of this model at early times depends or 
the value of a. 

(a) For a > 0, we have p gcg — > — oo when p gcg — > at 
z — Zb- Hence, there is a pressure singularity at z — z^, 
with a — > oo, while the scale factor a and the expan- 
sion rate a remain finite. This is the so-called "sudden 
singularity" previously discussed by Barrow |23(. In the 
classification scheme of Nojiri, Odintsov, and Tsujikawa 
|24| . this is a Type II singularity. The GCG density p, 
the equation of state parameter w, and the deceleration 
parameter q for this case are shown in Figs. 1-3, respec- 
tively, as solid curves. 

(b) For > a > — 1, both p gcg and p gcg go to zero 
at z = Zf,. Although w — > — 00 at z — Zb, there is no 
singularity in this case. Afterwards the GCG behaves as 
a growing phantom component and then asymptotically 
approaches the de-Sitter regime (w — > —1) at late times. 
The values of p, w, and q are shown in this case in Figs. 
1-3 as dashed curves. 

This model has also been examined in some detail in 
Ref. [HI. Their main point of emphasis is the interesting 
fact that this model allows for a phantom equation of 
state w < — 1, while avoiding a future singularity. We 
note further that for > a > — 1, this model is free of 
either past or future singularities, while still allowing for 
a phantom equation of state. 

Case 3) < A s < l,a < -1 (Transient GCG). In 
this case, one has a de-Sitter regime (w —> —1) at early 
times, while p gcg asymptotically approaches pressureless 
dust (w — > 0) at late times. The GCG density p, the 
GCG equation of state parameter w, and the decelera- 
tion parameter q for this case are shown in Figs. 1-3, 
respectively, as dot-dash curves. 

If the GCG serves as dark energy in this case, the accel- 
eration is a transient phenomenon. Models with transient 
acceleration are desirable from the point of view of string 
theory, as the existence of future horizons in an eternally- 
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FIG. 3: As Fig. 1, for the deceleration parameter q as a 
function of the scale factor a, where we have taken £l gC g0 = 0.7 
and f2 m0 = 0.3. 



accelerating universe leads to a well-known problem in 
constructing the S-matrix in such models |25ll26|. Hence, 
a fair amount of effort has gone into constructing mod- 
els in which the currently-observed acceleration is a tran- 
sient phenomenon (^IllElllS- This case for the GCG 
represents another such model. 

Case 4) A s > l,a < -1 (Late Phantom GCG) In 
this case p gcg starts as a constant and then grows and 
eventually hits a singularity in the future. The equa- 
tion of state resembles that of a cosmological constant 
(w ~ —1) at early times, and then becomes phantom- 
like in the future. The future singularity occurs at 
1 + z s = [A S /(A S - 1)]^+^. Note that this singular- 
ity occurs at a finite value of the scale factor, which is 
different from the standard phantom scenario in which 
the scale factor blows up simultaneously with the energy 
density of the dark energy. This singularity also differs 
from that in Case 2, as in this case, we have p — > 00 and 
p — * — 00 as z — > z s ; in the classification scheme of Ref. 
[24j . this is a Type III singularity. The GCG density 
p, the equation of state parameter w, and the deceler- 
ation parameter q for this case are shown in Figs. 1-3, 
respectively, as dotted curves. 

It is well known that, under very general assumptions, 
dark energy arising from a scalar field cannot evolve from 
w < -1 to 10 > —1 or vice- versa |3l|; this has been 
dubbed the "phantom divide" [12 • The GCG model also 
displays a phantom divide, characterized by whether or 
not A 8 > 1 or A s < 1. For all models with A s < 1, we 
have w > —1 at all times and for all values of a, while 
A s > 1 gives phantom behavior (w < —1) at all times 



FIG. 4: Confidence contours in fi m o — A s (left) and a — A s (right) parameter space, marginalizing over a and Q,mo respectively. 
Solid and dashed curves are the 68.3% and 95% confidence levels, respectively. 



and for all a. 

Note also that all of the GCG models display asymp- 
totic de Sitter behavior (w — > —1). For a < — 1, the de 
Sitter behavior occurs asymptotically in the past, while 
a > — 1 gives de Sitter behavior in the asymptotic future. 
These results are independent of A s , which simply deter- 
mines whether w approaches —1 from above or below. 

Many other barotropic models have been discussed in 
the literature, e.g., the Van der Waals model [3^ and 
the wet fluid model |34|. However, the barotropic model 
which most closely resembles the GCG models examined 
here is the model of Ref. [35| (see also |24j), which has 
the equation of state 

p = -p-Bp? (13) 

While this model is qualitatively different from the GCG 
models, it displays similar behavior in certain limits, par- 
ticularly with regard to singularities in the evolution. In 
particular, for B > and (3 > 1, the density p is an in- 
creasing function of the scale factor, so that the second 
term in equation l|13fl is dominant at late times. Thus, 
this model approaches the behavior of our Case 4 model 
at late times, with a similar singularity [24ll35| . 



III. CONSTRAINTS FROM SUPERNOVA DATA 

In this section, we will examine constraints on these 
various versions of the GCG from the supernova la data. 
In deriving these limits, we assume that the GCG acts 
purely as dark energy, and we take the dark matter to be 
a separate component. 

The observations of supernovae measure essentially the 
apparent magnitude m, which is related to the luminosity 
distance 4 by 



where 



D L (z) = ^d L (z) , 
c 



is the dimensionless luminosity distance and 

d L (z) = (1 + z)d M {z) , 
with c?m(z) being the comoving distance given by 

1 



Also, 



d M (z) = c 



M=M + 5 log! 



-dz' 



o H(z') 



c/H \ 
1 Mpc J 



25 



(15) 



(16) 



(17) 



(18) 



where M is the absolute magnitude. 

For our analysis, we consider the data set compiled by 
Riess et al. [2a|. The total data set contains the previ- 
ously published 230 data points from Tonry et al. [3^ |. 
along with the 23 points from Barris et al. 38j . But Riess 
et al. have discarded various points where the classifica- 
tion of the supernovae was not certain or the photometry 
was incomplete, increasing the reliablity of the sample. 
Ultimately the final set contains 143 points plus the 14 
points discovered recently using the Hubble Space Tele- 
scope (HST), and this set of 157 points is named the 
"gold" sample by Riess et al. 

The data points in these samples are given in terms of 
the distance modulus 

p obs (z) = m(z) - M(z) = 51ogd L + 25, (19) 
where d^ is measured in Mpc. The x 2 is calculated from 



m(z) = M + 5 log 10 D L (z) 



(14) 



X 



E 



Mobs(Zi) — Hth(Zi] H ,C a ) 



(20) 




FIG. 5: The 68.3% confidence contours in a — A s parameter 
space for different values of fi m o- Solid, dashed, dash-dot and 
dotted lines are for QmO = 0.35, 0.30, 0.25, and 0.20, respec- 
tively 



FIG. 6: Evolution of the matter density perturbation S as 
a function of the scale factor a (normalized to a = 1 at the 
present) for the standard GCG case with fl m o = 0.3. 



where present day Hubble parameter, Hq, is a nuisance 
parameter and c a are the model parameters. Marginaliz- 
ing our likelihood functions over the nuisance parameter 
Hq yields the confidence intervals in the c Q parameter 
space. 

For our purposes, we have three model parameters, 
^ m o, oi and A s . For our best fit analysis, we vary the 
parameters in the following ranges: fl m o from 0.2 to 0.4, 
a from —4 to 4 and A s from to 5. The best fit values 
for the parameters in this case are: £l m o = 0.39, a — 
-3.87 and A s = 4.99 together with x 2 min = 172.79. In 
Figure 4 we show the confidence contours in the O m o— A s , 
and a — A s parameter space by marginalizing over a and 
ft m o respectively. In both cases, A s = 1 (ACDM) is 
rejected at the 95% confidence level, while the data favor 
A s > 1 (wo < —1). Moreover, from the allowed a — A s 
space, one can see that a large portion of the allowed 
region corresponds to Case 4 (late phantom GCG). This 
is consistent with previous analyses that suggest that the 
supernova data slightly favor w < — I at present |3f5j . 
There is also a small region at the 95% confidence level 
corresponding to Case 2(b) (early phantom GCG without 
an initial singularity). 

To see how the model parameters depend on Q m o, 
we plot in Fig. 5 the 68.3% confidence contours in the 
a — A s space for different values of £l m o- This figure 
shows that the allowed parameters for the GCG depend 
very sensitively on the assumed value for O m o; for this 
set of choices for fi m o, any of the four cases is possi- 
ble. For fi m o = 0.2, the allowed region falls under the 
Case 1, where GCG behaves as a dust-like dark energy 
component, asymptotically approaching a cosmological 
constant. For Q m o = 0.25, the data allow Case I as well 



as Case 3 where the acceleration of the universe is only 
a transient phenomena as it asymptotically approaches 
dust-like behavior. For £l m o = 0.3 the data allow Case 
2 (both a and b) where the GCG behaves as phantom 
dark energy at early times but approaches a cosmologi- 
cal constant at late times, as well as Case 4 where the 
GCG evolves from cosmological-constant-type behaviour 
to phantom behaviour. For ft m0 — 0.35 only Case 4 is 
allowed. Our results show that O m > 0.3 is necessary in 
order to have a phantom-like equation of state. This con- 
clusion is consistent with the recent results of Jassal et 
al. |40| . although they used a different parametrization 
for the dark energy equation of state. 

These results show that the value of Cl m o is crucial in 
determining which types of GCG dark energy are consis- 
tent with the supernova data. The constraint on Q m p in 
ACDM models from WMAP is O m0 = 0.29 ± 0.07 [4l|, 
which is also consistent with recent results from SDSS ob- 
servations, while the more recent 2dFGRS analysis gives 
ft m0 = 0.237 ± 0.020 H2. Of course, when w is allowed 
to vary from w = — 1, these limits become functions of w, 
but it is clear that the current observational constraints 
on O m o are insufficient to rule out any of our four possible 
GCG models using supernova data alone. 

IV. GROWTH OF LINEAR DENSITY 
PERTURBATIONS 

In this section, we study the growth of density pertur- 
bations for the mixture of a matter fluid and a GCG dark 
energy fluid in the linear regime on subhorizon scales. In 
performing this calculation, it is necessary to assume a 
particular clustering behavior for the dark energy. How- 
ever, this behavior will depend on the physical model 
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FIG. 7: Evolution of the matter density perturbation 5 as 
a function of the scale factor a (normalized to a = 1 at the 
present) for the transient GCG case with Q m o = 0.3. 
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FIG. 8: Evolution of the matter density perturbation S as 
a function of the scale factor a (normalized to a = 1 at the 
present) for the late phantom GCG case with Q m o = 0.3. 



that gives rise to the Chaplygin gas equation of state. 
For example, if the Chaplygin gas is taken to be a perfect 
fluid satisfying equation JJjJ, then the GCG component 
will cluster gravitationally with a sound speed given by 
c 2 s = —wa 7]. On the other hand, it is also possible to 
generate minimally-coupled scalar field models with the 
equation of state given by equation |JBJ 4|, |2l| (see also 
Ref. [13). Such models always have c s = 1 on subhori- 
zon scales and therefore do not cluster on small scales. 
The evolution of density perturbations for a GCG dark 
energy component that can cluster gravitationally was 
examined in Ref. , while the case of a smooth, unclus- 
tered GCG dark energy component was examined in Ref. 
[iflj . We take the latter approach here. In this case, the 
only effect of the GCG evolution is to alter the growth of 
dark matter perturbations through the effect of the GCG 
energy density on the expansion of the universe. How- 
ever, it is important to remember that the results will be 
very different for the case in which the GCG component 
is allowed to cluster. 

Assuming the GCG to be a smooth component, the 
growth equation for the linear matter density perturba- 
tion, 8, is given by 

S" + (2 + — ^)S' + 3 Cl S = 0, (21) 

where "prime" denotes the derivative with respect to 
ln(a), "dot" denotes the derivative with respect to t, and 
H is the Hubble parameter for the background expansion 
given in equation i|10[l . In equation (|21(l . S is the linear 
matter density contrast, S = 5p m /p m , and c\ is given by 



2 Q m0 + {l gcg0 [l + A a ( a 3(°+D - l)]i/(i+«) ■ 



One can easily check that for A s = 1, the equation re- 
duces to that for the ACDM model. We have integrated 
equation l|21|) numerically from a = 10~ 3 to a = 1 (taken 
to be the present). The initial conditions are choosen 
such that at a = 10 -3 , the standard solution S ~ a for 
Einstein-deSitter universe is reached. Also we have as- 
sumed the matter density parameter Q m o = 0.3 through- 
out. We have studied the solutions for the standard GCG 
(Case 1), transient GCG (Case 3) and late phantom GCG 
(Case 4) models. 

Figs. 6, 7, and 8 show the behaviour of 5 as a function 
of the scale factor for these three cases. The standard 
GCG case has been studied previously by Multamaki et 
al. |43j , who showed that for parameters slightly deviat- 
ing from the ACDM universe (A s = 1), 6 deviates grossly 
from the standard ACDM result, making this model hard 
to reconcile with the present observational data. The be- 
haviour of S in Fig. 6 agrees with this result. 

On the other hand, one can see from Figs. 7 and 8 (the 
transient and late phantom cases respectively), with val- 
ues of a and A s deviating significantly from the ACDM 
case, that the behaviour of 6 is practically indistinguish- 
able from the ACDM case. This is an interesting result; it 
shows that models with GCG dark energy can have quite 
different equations of state from ACDM at either early 
or late times, yet still give similar results for the growth 
of linear density perturbations. The reason for these re- 
sults is clear from the behavior of p for these models (Fig. 
1). Both the transient and late phantom GCG models, 
like the cosmological constant, contribute negligibly to 
the density of the universe at early times; this density 
is dominated by the matter component. At low rcdshift, 
the GCG in both models begins to dominate the expan- 
sion (just as the cosmological constant does in ACDM), 
with the GCG density decreasing in time (for the tran- 
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FIG. 9: Comparison of the evolution of the matter density 
perturbation S for the transient and standard GCG cases with 
SimO = 0.3. 

sient model) or increasing in time (for the late phantom 
model). However, these deviations from the behavior of 
the cosmological constant occur over a very short range 
in redshift, and by forcing ilo for the dark energy to be 
the same in all three cases, the results for the evolution 
of density perturbations are almost exactly the same. 

For the standard GCG, in contrast, the dust-like be- 
havior results in a significant contribution to the density 
of the universe as long as dark matter is the dominant 
component; if the GCG is assumed not to cluster, the re- 
sult is a significant decrease in the perturbation growth. 
This is illustrated more clearly in Fig. 9, where we show 
5 as function of the scale factor for both the transient 
and standard GCG models with the same A s (A s = 0.7), 
along with the ACDM model. 

Since the observations related to density perturbations 
are consistent with the ACDM model, our results suggest 
that, unlike the standard GCG model, both the transient 
and late phantom GCG models are consistent with the 
linear growth of density perturbations inferred from ob- 
servations. We have not shown the results for the early 
phantom model, but for the case where there is no early 
singularity, the results are also nearly identical to the 
ACDM model. Again, we emphasize that these results 
assume a non-clustering GCG. For the case where the 
GCG clusters as a perfect fluid, the growth of density 
perturbations would be quite different. 

V. CONCLUSIONS 

Our exploration of the the full A s — a parameter space 
for the generalized Chaplygin gas yields three additional 



models beyond the standard GCG. Although these can- 
not serve as unified models for dark matter and dark 
energy, they can have interesting consequences when 
treated as models for dark energy alone. 

The early phantom model (Case 2) has the interesting 
property of serving as phantom dark energy (w < —I) 
without a late big rip singularity, and for some choices 
of parameters it is also free of an initial singularity. The 
transient CGC model (Case 3) provides a mechanism for 
accelerated behavior at the present, but it asymptoti- 
cally approaches dust-like behavior at late times (i.e., 
its time evolution is exactly opposite to the standard 
GCG model). Thus, it provides a mechanism to allow for 
present-day acceleration without a future horizon. The 
late phantom GCG model gives w < — 1 with w decreas- 
ing with time, and it results in a future singularity at a 
finite value of the scale factor. 

All of these models can be made consistent with 
the type la supernovae observations, for an appropriate 
choice of Q m o; the question of which models are allowed 
is extremely sensitive to the value of f2 m o- If the GCG 
is assumed not to cluster, then all of these models (with 
the exception of the subset of early phantom models with 
an initial singularity) are also consistent with the growth 
of linear density perturbations. Of course, if the GCG is 
assumed to cluster, then these results will be significantly 
altered. 

Note that all three of our "new" GCG models have a 
de Sitter phase, and all three models can be tuned ar- 
bitrarily close to a cosmological constant at the present, 
either by pushing the phantom-like behavior arbitrarily 
far into the past (early phantom model), or by pushing 
the dust-like or phantom-like behavior arbitrarily far into 
the future (transient GCG and late phantom model, re- 
spectively). These limits may seem uninteresting, as they 
reduce to the ACDM model over all observable ranges, 
but the one exception is the transient GCG model. A 
dust-like phase for this model, even in the far future, 
eliminates the problem of future horizons. Thus, the 
transient GCG model can be made arbitrarily similar to 
the ACDM model but at the same time can resolve the 
possible conflict between the accelerating universe and 
string theory. 
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